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We point out that by modifying the setup of a recent experiment that generates a Dirac String (4), one can 
create a quasi 2D spinor Bose condensate on a cylindrical surface with a synthetic magnetic field pointing 
radially outward from the cylindrical surface. The synthetic magnetic field takes the form of the Landua gauge. 
It is generated by the Berry’s phase of a spin texture, frozen by an external quadrupolar magnetic field. Unlike 
in the planar case, there are two types of vortices (called A and B) with the same vorticity. The ground state for 
5 < S < 9 consists of a row of alternating AB vortices lying at the equatorial circle of the cylinder. For higher 
values of S, the A and B vortices split into two rows and are displaced from each other along the cylindrical 
axis z. The fact that many properties of a BEC are altered in a cylindrical surface implies many rich phenomena 
will emerge for ground states in curved surfaces. 


In the study of quantum matter, one usually deals with Eu¬ 
clidean space. Spaces with non-zero curvatures are seldom 
encountered. Yet in theoretical studies, torus and spherical 
surfaces are often used for computational convenience, or to 
demonstrate hidden topological properties of the system. In 
recent years, it is found that curvature effects can mimic those 
of gauge fields (as in the case of graphene (TJ), and can lead to 
topological responses such as Hall viscosity (2j. The fact that 
spatial curvature can help uncover new features of quantum 
matter makes it desirable to create manifolds of controllable 
curvature, and to develop capability to add in synthetic gauge 
fields. 

The purpose of this paper is to point out ways to create 
quantum gases in curved surfaces and to manufacture syn¬ 
thetic gauge fields in such systems. As a first step, we shall 
consider the simple geometry of an annulus and the case of 
spinor condensate of bosons with large spin. We consider 
Bose-Einstein condensates (BEC) because its quantum prop¬ 
erties are often magnified considerably by Bose statistics. 
Large spins particles are considered because the Berry phase 
that determines the strength of synthetic gauge fields for these 
particles is proportional to their spin. As we shall see, the sys¬ 
tem we construct amounts to a system of ’’charged” bosons in 
a magnetic field written in Landau gauge. The effective mag¬ 
netic field pointing radially outward, normal to the cylindrical 
surface. This gauge field is caused by a spin texture generated 
by a (real) quadrupolar magnetic field. We shall also point out 
that the techniques we apply in our construction have all been 
achieved experimentally; therefore, creating BEC on cylindri¬ 
cal surface with a Landau gauge is feasible. 

As it turns out, for typical parameters of atomic gases, a 
ground state in the lowest Landau level (“mean field quantum 
Hall regime”) can only be achieved with a thousand of par¬ 
ticles or less. A Bose gas with with 10 5 bosons will form a 
condensate occupying many Landau levels. Lor a sufficiently 
strong synthetic magnetic field, the ground state will contain 
vortices El . However, due to the topological constraint of the 
annulus and the curvature of the cylinder, the flow pattern of 
a single vortex as well as the vortex lattice are very different 
from those found in planar geometry. Lor a BEC in a plane, 


isolated vortices have cylindrical symmetry and the vortex ar¬ 
ray of a rotating BEC is a hexagonal lattice. Lor BECs in 
an annulus, however, there are two kinds of vortices with the 
same vorticity but different flow pattern. In a large effective 
magnetic field, the vortices will first align in a single row at 
the center of the annulus at zj = 0 in the azimuthal direction ip. 
After certain proliferation, the vortices split into two lines at 
Zj = ±Z. The velocity profile resembles that of two counter- 
circulating superfluid rings separated in the z direction. In the 
time of flight experiments, a row of n vortices will lead to an 
ft-fold symmetry in the density profile in the azimuthal direc¬ 
tion. The presence of many quantum phenomena have taken 
new forms in curved space even for the simple case of Bose 
condensates suggest a great deal more new phenomena lay in 
store for more complex quantum systems. 

I. Experimental setup and realization of the Landau gauge 
in a cylinder : Our setup is a modification of the recent ex¬ 
periment by Daivd Hall’s group to generate a Dirac string in 
a BEC by inserting in it a quadrupolar magnetic field|4|. Our 
system is shown in Lig. [I] A confining potential of the form 
an annulus is constructed by piercing through an attractive po¬ 
tential (produced by a red-detuned laser) with a repulsive core 
(produced by a blue-detuned laser). This will create a confin¬ 
ing well in the radial direction with a minimum at radius R. 
A harmonic potential V(z) = Moj^z 2 / 2 is imposed along the 
z-direction, which is far weaker than the harmonic trap in the 
radial direction in the neighborhood of R. When a quantum 
gas is loaded onto this trap, it will form a cylindrical layer 
with radius R and thickness cr. By increasing R and reducing 
the frequency cj Z9 we shall have R » cr. The quantum gas 
then becomes a curved quasi-2D system. 

Next, we insert a quadrupole magnetic B at the center of the 
cylinder, 

B(x) = B 0 (xx + yy- 2 zz) 

= #oKsin 6 cos ipx + sin 6 sin ipy - 2 cos 6z :). (1) 

The direction of B on the cylindrical surface is shown in Lig. 
[T] This quadrupole field is introduced to generate an effective 
gauge field for the bosons that are spin aligned with the lo¬ 
cal field. The insertion of a quadrupolar magnetic field into 
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FIG. 1: Schematic plot of the system. 


a 3D BEC has been recently performed in the experiment in 
Ref.El. Our configuration is basically a modification of the 
set up in Ref. (4) there by piercing through their harmonic po¬ 
tential with a repulsive blue-detuned laser. For a sufficiently 
large quadrupole field, the spins of the bosons will align with 
the local magnetic field. Denoting direction of the spin as 
/ = cos a\z + sin a(cos/3i: + sin/3y), Eq.{l]) implies 


a - ip. 


cos/3 = 


2 z 

VR 2 + 4 z 2 


( 2 ) 


The condensate wave function of bosons with spin S is then 
0^(x) = £ fl (x)0(x), where a is the spin index, £ a (x) is a 
normalized vector aligned with the local magnetic field, i.e. 
B(x) • S a b£b( x ) = S<Ta( x )- With the spin direction frozen by B, 
the kinetic energy becomes 


|V^| 2 = 




+ \vtf\<t >\ 2 + (pvtfm 2 - 


(3) 


The field 0 experiences an effective gauge field A(x) = 
-/^(x)Vf(x). The energy functional of 0 (with spin indices 
suppressed), 


-I- 


e w= i 


(M-v( Z ))w 2 + ^m 4 


(4) 


reduces to a functional of 0, and the last two terms in Eq.Q 
(which we shall see is only a function of z) can be absorbed 
into potential V(z). In the limit where the radius R of the quan¬ 
tum gas is much larger than its thickness <x, the system can be 
regarded as a quasi-2D system, and the field 0 is a function of 
z and the azimuthal angle ip. It is then straightforward to show 
that 


V0 = zd z (p + -d ^0, 

= -S(Va) cos/3 = 


0 2 zS 


0 2 zS 


(5) 

( 6 ) 


R VR 2 + 4z 2 R R 
V X A = (r ~ 1 d, f A- - d z A, f ) f = -(2S/R 2 )r, (7) 


for z < R/ 2, where r is the radial unit vector in cylindrical 
coordinate. This effective magnetic field corresponds to 2 S 
unit of phase winding in the planar case. Thus, for Bose gases 
with sufficiently large spin, the ground state is expected to 
contain vortices. 

To simplify notations, we measure length in units of R , 
so that z/R —» z, and z is now dimensionless. We further 
introduce the dimensionless variables a = McJj z R 2 lfi,p = 
p/(% 2 /2MR 2 ),g = g/(fi 2 /2MR 2 ), then the energy functional 
becomes 

ew = 2 

S = \d z (f>\ 2 + | [~id v + 2 S z) <j >| -{fi- a 2 z 2 ) \</>\ 2 + ||^| 4 . (8) 

The first two terms represent a system in uniform magnetic 
field with a Landau gauge. 

II. Isolated vortex and vortex array patterns on a cylinder: 
To find the vortex pattern for the cylindrical 2D condensate, 
we seek for a variation solution in the Thomas-Fermi (TF) 
limit for the functional Eq.([8|. We use the variational anzatz 

<P = yfnfe 1 ® ( 9 ) 


Here n(yp, z) is the density profile given by the TF approxima¬ 
tion in the absence of vortices. f(<p,z) punches hole in the 
vortex locations (<pi, Zi ) and modifies the density profile, 


= [4 Lanli 


V(y - <ft) 2 + (z - Zi ) 2 


( 10 ) 


where £ is the core size, also written in units of R. The phase 
0(<£>, z) gives a 2n phase winding around each vortex. This 
form has been shown to match well with experiment in the 
case of rotating gases 01. 

The periodicity of the wave function (or e 1 ®) in cp, however, 
strongly constraints the form of 0. From the conformal map¬ 
ping between a punctured plane and a cylinder, one finds two 
distinct phase factors e l ®j that have the same +2 n phase wind¬ 
ing about the local (ypj, zj ), 


exp(/0p = |^|, W ± = e ±iu - e ±m > (11) 


where u = cp + iz , and uj - (pj + izj is the vortex location. The 
superfluid velocity of each vortex is ^ V0* where 


d v &j = t-t - 1 — 

2 [cosh(z - Zj) - cos (ip - (pj 


= 


_ sin (tp-tpj) _ 

2 [cosh(z - zj) - cos (<p - (pj)) 


( 12 ) 

(13) 


One can easily identify that V0t is related to V0T by a n 
rotation with respect to the vortex core. For |z - Zj\ » 0, they 
approach 


V0t 


-i-^z»Zj 

0, z « Zj ’ ' 


0, z » Zj 
1 • <p, z «: Zj 


( 14 ) 
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So, as a function of z, <9^0 j (or <9^07) approaches to a con¬ 
stant far above (or far below) the vortex core, and vanishes 
on the other side, as shown in Fig. 2a and Fig. 2b Because 
of this feature, we call 0t and 07 the A- and B-vortex as 
their superfluid velocities are mostly non-vanishing ’’above” 
and ’’below” the vortex core respectively. Note that despite 
their very different velocity patterns, both of them have the 
same +2n circulation around the vortex core. In the case we 
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(c) 6 alternating vortices aligned along z = 0. 




FIG. 2: Velocity profiles for two types of vortices and the 
vortex array. 


have A-vortices on a set of Q point {uj,j = 1 ,..Q] and B- 
vortices on another set of Q points, {uy, j = Q + 1, ..2 Q], the 
phase factor is 


Q 2 Q 

exp(/0) = W/\W\, W = fk n Wj,. (15) 

i= 1 7 = 0+1 

We have minimized the energy Eq.([8]) with the variational 
wave function Eq.Q, where the vortices at (ipj, zj) can be ei¬ 
ther A or B type. We shall first present our results and describe 
the details of our calculation at the end. The variational result 
is shown in Fig. [3] First, since the system has reflection sym¬ 
metry in z-direction, A and B vortices must appear in pairs at 
appropriate location to respect this symmetry. Our variational 
calculation shows that for S < 4, the strength of the gauge 
field is not strong enough to generate vortices in the ground 
state. For S =5,6,7,8, there are 4,6, 8,10 vortices respec¬ 
tively lying on the circle at z = 0 (i.e. zj - 0). These rows 
of vortices are all in the alternating pattern A-B-A-B-.. . For 
S > 9, the vortices split into two rows above and blow z = 0, 
with the A-vortices shifted up and B-vortices shifted down. 

In our calculations, we have considered a system of 10 5 
bosons, with a chemical potential such that the extent of the 
condensate along z is of the order but shorter than R to stay so 
that A has the form of the Landau gauge, Eq.([6]). The size of 
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FIG. 3: Number of vortices versus spin in the ground state. 

For S < 4, the ground state contains no vortex. Within 
S = 5,6,7,8, the vortices aligned in one row at Zi = 0 with 
the pattern A-B-A-B. For S > 9, the vortices array splits into 
two rows centering at Zi = ±Z, with all A vortices aligned in 
one row and all B vortices in another. 


the vortex core is chosen to be the coherence length of uniform 
condensate with the chemical potential at z = 0. 

We conclude this section by examining the phase function 
and velocity field of the alternating vortex row in greater de¬ 
tail. Let us consider the case of S = 6 (corresponding to 168 Er) 
where the ground state has 6 vortices of alternating A and B 
type equally spaced on the circle at z = 0, say, at (p = 
where n - 0,1, ..5. Defining w = e lu = e l(p ~ z , and a = e 2n ' 6 , 
the phase function W in Eq.(p~5]) is then 

W= n [o - a ln )(w~ l ~ a- (2n+l) )] (16) 

11=0,1,2 

= w 3 - w" 3 (17) 

The wave function is simply 

- e 3i<p-3z _ e -3i<p+3z 

<P(<p,z) ~ VT r 7/(y,z) |g 3 ^ _ e _ 3lV+3z| , (18) 

Eq.(p~8]) shows that for z > 0 (z < 0), <f>((p,z) quickly ap¬ 
proaches e~ 3l(p ( e 3l(p ). The system then consists of two counter 
circulating superflow above and below z = 0, as shown in 
Fig.2(c). 

The reason that the vortex pattern is so different from the 
planar case is a consequence of the confined geometry. As one 
increases the strength of the gauge field by, say, increasing S , 
more vortices will come. At some point, it will be too costly 
to put all the vortices in a single row. Our calculation shows 
that the vortices will separate into two rows, as shown in Fig.3 
for the cases S > 9. Mathematically, if one squeeze in a large 
number of vortices by choosing a very large value of S , the 
vortices will then organize in a lattice on a cylinder as in the 
planar case. The separation of vortices in two rows as shown 
in Fig.3 can be viewed as a tendency towards the lattice limit. 

III. Signature of the vortex array: The presence of these al¬ 
ternating vortex array can easily be detected in time of flight 
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FIG. 4: Density profile at xq /R = y/htJm/R = 7 on the x-y 
plane for z = 0. Length is measured in units of R. 


experiments. As we have discussed, the AB vortex row con¬ 
sist of 2 n vortices will generate two counter phase superflow 
e~ in(p and e +in<p for z > 0 and z < 0. The system can be approx¬ 
imated by two rings of condensates with opposite circulation, 
with one ring sitting above the other along z with a separation 
or the order of their radius R. In the time of flight experiment, 
these two rings will produce an interference at the z = 0 plane 
of the form of e~ incp + e +irup and therefore exhibit a density pat¬ 
tern with 2^-fold symmetry. This effect is in fact found in an 
explicit calculation of the ballistic expansion of the vortex row 
condensate in Eq. ©• The time evolution of the condensate 
is given by 


0(x,o= f d Z r 

J -1 J-7T 


d(pU((p , z, t; (p\ (19) 


where U((p 9 z,t',<p',z?) is the Green’s function for free 
particle propagation at large distance and at long 
times in cylindrical coordinates, U((p,z, t; (p',z f ) ~ 

exp [-z(Rrcos(<£ - ip') + zz')/*o]> and xq = ^JWfm. The 
density pattern at the equatorial plane z = 0 at long times is 
shown in Fig.4. 

IV Details of the variational calculation: Since both n , / in 
Eq.Q are real and \h\ = 1, we can rewrite Eq.([ 8 ]) as 

S — 8\ + 82 + 83, 


fii = [(d : 0) 2 + (c\-0 + 2Szf\nf 


f 


1 


gif 


62 = _ «/V / - - J - - -) - V 


[<v V «) 2 


0 u- 


V( z))n + | n 2 


r 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


We recognize 83 has the form of usual Gross-Pitaevskii en¬ 
ergy functional, and determines the density profile in the ab¬ 
sence of vortices. Applying TF approximation to Equation 
(23) and vary with respect to n , we obtain 

n - V(z) Moj 2 R 2 


n(z) = 


g 


2 g 




(24) 


with the constraint n > 0. Here /3 = 




, and the 


condensate spreads within Taking the condensate to 

be homogeneous along radial direction, the total number of 
particles given by this density profile is then 



a s /cr’ 


(25) 


where a s is the scattering length related to g by g = 71 M as . 

Further, in the absence of vortices, / = 1 everywhere, and 
the first two terms in 82 vanishes. That means 82 is the en¬ 
ergy increase due to the existence of vortices. To analyze it, 
we first notice that the core size is given by the coherent length 
f ~ , n 2 Taking the value at z = 0, we have £0 ~ A • 

y 2 Mn(z)g K a 

Note a = R 2 /d 1 , where d = is the natural length of the 
harmonic trap along z, which we take to be an order of mag¬ 
nitude smaller than R. Then f ~ 10 2 « 1. That means 
the first two terms in 82 is only non-zero close to each vor¬ 
tex core. So we can make approximations using the identity 
£* dxfi(x)f 2 (x) = (Zifi(xi))fZ,dxf 2 (x), where /i (x) is a 
smooth function, and fz(x) is only non-zero within a small 
region around xf s. Then the integration can be done analyt- 
ically: f d 2 x S 2 = + (2 + a 2 f-)C ~ z 2 ), where 

C = |(41n2 - 1) « 0.928125 is the integration result, and Zi 
is the z coordinates of vortices, and it is confined to \zt\ < 

As expected, the vortex core energy would like to repel vor¬ 
tices out of the condensate by increasing \zf The constant 
term plays no role and will be omitted. Also, we take af - 1 
according to the above analysis. 

In summary, measured in units of Mo ^ R , the energy func¬ 
tional is reduced to 


e [, A] = 3cyyi-z 2 )+ 

i 

X [(cLQ) 2 + 


rn rP 

I dip f dz 
J-7T J-P 

(d v ® + 2 Sz) 2 ](J3 2 -z 2 )f 2 . (26) 


For fixed total particle number N, reducing /3 corresponds to 
increasing the strength of the trap or weakening the interac¬ 
tion. We will take = 1 in the following. The vortices 
should induce the superfluid velocity profile that serves to can¬ 
cel (2 S z) with <9^0 in the above energy functional, and there¬ 
fore lower the energy. 

Concluding Remarks: The emergence of two kinds of vor¬ 
tices with identical vorticity in a cylindrical manifold is a new 
feature of Bose condensates in a cylinder. It is a consequence 
of the topological constraint on the single value-ness of the 
wave function (i.e. that forces the spatial dependence to be 
expressed in terms of e in(p ) and should persists even when the 
manifold is deformed. Although we focus on a particular as¬ 
pect of the quantum gas in curved surfaces, there are a lot 
more to explore especially for systems with more complexity. 
Realization of quantum gases in curved surfaces will surely 
open an exciting direction for cold atom research. 
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